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ABSTRACT 


Through  an  analysis  of  a  more  dimensional  extension  of  the  bulk  mixed 
layer  model  of  Kraus  and  Turner  (1967)  it  is  shown  how  even  under  spatially 
uniform  atmospheric  conditions  an  initially  smooth  horizontal  temperature 
gradient  in  the  surface  mixed  layer  can  still  develop  into  a  front.  Essential 
for  it  is  the  notion  that  the  net  downward  heat  flux  at  the  air-water 
interface  is  related  to  the  difference  between  the  mixed  layer  temperature 
T  and  an  apparent  atmospheric  temperature  Tft  (Haney,  1971)  so  that  the 
initial  horizontal  gradient  in  T  corresponds  to  a  horizontal  variation  in 
the  surface  buoyancy  flux.  As  a  result  also  the  mixed  layer  depth  differs 
from  place  to  place.  Depending  on  the  direction  of  the  wind  driven  transport 
this  produces  either  a  steepening  or  a  flattening  of  the  initial  temperature 
profile.  A  lower  bound  condition  for  the  initial  horizontal  advective  heat 
flux  is  derived  in  terms  of  the  initial  surface  heat  flux,  the  stirring  by  the 
wind  and  the  time  rate  of  change  of  the  apparent  atmospheric  temperature.  If 
that  condition  is  satisfied  a  front  develops.  If  not,  then  frontogenesis  is 
prevented  by  the  damping  effect  of  the  local  atmospheric  heating.  It  is  shown 
that  the  condition  can  be  satisfied  on  the  scales  of  lakes  (or  small  seas)  and 
in  near  equatorial  oceanic  regions.  Mathematically  the  problem  boils  down  to 
a  first  order  quasilinear  hyperbolic  partial  differential  equation  that  is 
solved  exactly  by  the  method  of  characteristics.  , 


AMS  (MOS)  Subject  Classifications!  35L67,  76L05,  86A05 

Key  Wordst  frontogenesis,  mixed  layer  dynamics,  air-sea  interaction, 
quasilinear  heat  equation 

Work  Unit  Number  2  (Physical  Mathematics) 


* 

Cooperative  Institute  for  Marine  and  Atmospheric  Studies,  4600  Rickenbacker 
Causeway,  Miami,  FL  33149.  (On  leave  of  absence  from  Rijkswaterstaat,  Data 
Processing  Division,  Ni jverheidsstraat  1,  2288  BB  Rijswijk  (z.h.).  The 
Netherlands . ) 


This  work  was  supported  by  the  National  Oceanic  and  Atmospheric  Administration 
-  Environmental  Research  Laboratory  under  Contract  No.  NA79  RAD0001 0 >  and 
partially  sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C- 
0041,  and  the  Netherlands  Ministry  of  Transportation  and  Public  Works. 


cor* 

msricreP  , 


FRONTOGENESIS  IN  AN  ADVECTIVE  MIXED  LAYER  MODEL 


Will  p.  M.  de  Ruijter 


NT IS  GRAAI 
r'nc  tas  □ 

U  iiiuouncsd  □ 

Justification--. _ 


|  Distribution/ 

Availability  Codas 
Avail  and/or 
jDlst  Spatial 


1.  Introduction  J!  \ 

Until  rather  recently  bulk  models  developed  for  the  upper  mixed  layer  of  the  sea  were 
mainly  one-dimensional  (see  e.  g.  the  review  by  Hiiler  and  Kraus,  1977).  This  was 
motivated,  among  other  things,  by  studies  like  the  one  by  Gill  and  Miller  (1973)  who  used 
rough  scaling  arguments  to  show  that  advection  plays  only  a  minor  role  in  the  evolution  of 
the  mixed  layer  in  most  oceanic  regions  not  too  close  to  the  equator.  Wear  the  equator 
wind-driven  transports  in  the  upper  layer  can  be  considerably  larger  so  that  horixontal 
advection  can  play  a  more  prominent  role  in  the  mixed  layer  dynamics  of  those  regions. 

Also  in  stratified  lakes  and  reservoirs,  with  their  completely  different  scales,  river 
throughflows,  selective  withdrawal  and  intrusions,  possible  boundary  and  bottom  effects 
etc.  (see  the  review  by  Imberger  and  Hamblin,  1982)  advection  can  play  a  dominant  role. 

One  of  the  most  dramatic  effects  of  horitontal  advection  can  be  the  formation  of 
fronts  in  the  upper  layer  and  it  is  mainly  in  this  context  that  the  one-dimensional  bulk 
mixed  layer  models  were  recently  extended  to  the  more  dimensional  case.  Examples  are  the 
steady  state  model  of  Welander  (1981)  and  studies  by  De  Ssoeke  (1980)  and  Cushman-Roisln 
(1981).  In  all  these  models  the  physical  process  that  is  responsible  for  the  generation  of 
fronts  is  the  horizontal  variability  of  the  wind  stress  such  that  it  produces  a  convergent 
Ekman  transport. 

In  this  study  we  will  isolate  a  completely  different  physical  mechanism  that  can 
produce  frontogenesls.  It  does  not  depend  on  a  horizontal  gradient  in  either  wind  stress 
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or  atmospheric  temperature  at  the  surface;  both  shall  be  spatially  uniform  throughout  this 
papar.  An  assantlal  ingredient  for  It  la  that  the  downward  heat  flux  at  the  surface  is  not 
merely  a  prescribed  function  but  related  to  the  difference  between  the  mixed  layer 
temperature  and  an  apparent  atmospheric  equilibrium  temperature  TA  that  includes  the 
effects  of  evaporation  and  solar  radiation  (Haney,  1971).  An  initially  smooth,  horizontal 
gradient  of  the  mixed  layer  temperature  thus  produces  a  variability  in  the  surface  buoyancy 
flux.  Associated  with  that  also  the  mixed  layer  depth  shows  a  horizontal  variation,  with 
a  uniform  wind-driven  transport  this  can  under  certain  conditions  (derived  in  Section  3) 
lead  to  a  steepening  of  the  temperature  profile  with  eventually,  in  the  absence  of  a 
horizontally  diffusive  process,  the  formation  of  a  discontinuous  step  in  the  temperature 
distribution. 

In  Section  2  we  present  the  basic  mathesutical  formulation  of  the  model  which  is  in 
essence  the  sure  dimensional  extension  of  the  classical  nodal  by  Kraus  and  Turner  (1967). 

In  Section  3  we  then  show  that  for  the  non-entraining  application  the  process  can  be 
described  by  a  first  order  quasilinear  hyperbolic  partial  differential  equation  that  can  be 
solved  analytically  by  classical  method*.  Inspiration  for  its  analysis  (section  4)  can  be 
based  on  the  study  of  related  problems  in  completely  different  fields  of  application 
ranging  from  shock  waves  in  gas  dynamics  (Courant  and  Friedrichs,  1948)  to  the  evolution  of 
jams  in  traffic  flow  (Light hi 11  and  Whit ham,  1955;  Richards,  1956).  Many  of  those  specific 
applications  are  described  in  Hhitham  (1974). 


2.  Formulation  of  th«  modal 


In  this  section  the  mathematical  formulation  of  the  advective  mixed  layer  model  shall 
be  presented.  The  starting  point  is  the  heat  conservation  equation,  in  which  horisontal 
diffusion  is  omitted: 
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T  is  the  bulk  temperature;  u,  v  and  w  are  the  velocity  components  along  the  x,  y 
and  a  axes,  respectively;  x  is  positive  eastward,  the  positive  y  axis  points  northward 
and  a  is  positive  upward  from  the  watarsurface  t  *  0|  w'T'  is  the  vertical  component 
of  the  diffusion  flux  density;  primes  denote  deviations  from  the  mean,  overbars  denote 
(ensemble)  averages.  PQ  is  a  reference  density  and  c  the  specific  heat.  I  stands  for 
the  penetrating  component  of  solar  radiation;  I  ■  I^e^*,  with  Iq(<  0)  the  surface  flux 
of  penetrating  solar  radiation  and  Y  an  extinction  coefficient,  ranging  from  0.3m-1  in 
polluted  Dutch  lakes  to  0.03m"1  in  certain  areas  in  the  ocean. 

With  approximately  homogeneous  turbulence  in  horisontal  planes  and  neglecting  the 
production  of  turbulence  kinetic  energy  by  the  mean  flow  shear  the  steady  state  equation 
for  turbulent  mechanical  energy  is  (cf.  Phillips,  1977) 


f-  (w*  <£-  ♦  i  (u'2  ♦  v'2  ♦  w’2))}  -  ogw'T'  -  e  .  (2) 

**  p0  * 

It  expresses  a  balance  between  the  divergenca  of  the  turbulent  vertical  energy  flux,  the 
rate  of  working  by  the  buoyancy  force  and  the  dissipation  ( e) .  u  is  the  thermal 
expansion  coefficient  and  p'  denotes  the  pressure  fluctuations. 

The  assumption  that  the  mean  temperature  is  uniform  over  the  depth  of  the  mixed  layer 
makes  it  easy  to  integrate  the  governing  equations  in  the  vertical  over  the  mixed  layer. 
For  the  one-dimensional  case  this  was  for  Instance  worked  out  and  discussed  in  detail  by 
Miller  and  Kraus  (1977).  A  more  dimensional  extension  was  systematically  derived  by  De 
Szoeke  (1980).  Except  for  the  addition  of  advective  terms  in  the  heat  transport  equation 
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an  Important  dlf faranca  with  tha  one-dimensional  models  la  that  now  tha  notion  of  tha 
intarfaca  batwaan  tha  mixed  layar  and  tha  intarior  can  also  ba  controllad  by  up-  or 
downwelling  (vortical  advection).  Tha  resulting  vertically  integrated  fora  of  tha  haat 
tranaport  aquation,  (1),  iai 


h  H  +  U  5  +  V  f  +  tT  -  Ti>*. 


-w*T* (0 )  -  ^ 


whara  h  ia  tha  aixad  layar  depth, 


(U,V)  -  /  (u,v)da,  tha  horiaontal  tranaport 
-h 


(3) 


w'T'(O)  tha  aurfaca  haat  flux  and 

Tj  tha  taaparatura  iaaodiataly  below  tha  aixad  layar. 

Tha  antrainaent  velocity  wa  ia  given  by  tha  relation 

Sh  ,  3U  ,  3v  ... 

%  “  at  ax  ay  l4> 

axpreaaing  that  in  the  caaa  of  antrainaumt,  i.a.  wa  >  0,  the  rata  of  change  of  the  aixad 
layer  depth  ia  a  balance  batwaan  antrainaent  and  tha  divergence  of  the  tranaport  in  tha 
aixad  layar.  In  the  detraining  caaa  (wa  <  0)  tha  aixad  layar  la  decoupled  from  tha 
interior  aa  there  ia  no  temperature  jump  at  tha  baae  of  tha  mixed  layer. 

Vertical  integration  of  (2),  neglecting  the  flux  of  mechanical  energy  at  the  base  of 
the  mixed  layer,  givoa  the  mechanical  energy  balancet 
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whara  u*  =  I i/pQ |  *  ia  tha  friction  velocity  (with  t  tha  surface  wind  atrass)  and  a 
a  "mixing  efficiency"  parameter.  Tha  mechanical  energy  balance  states  that  that  portion  of 
tha  kinetic  energy  input  by  tha  wind  that  is  not  disaipated  is  used  to  increase  tha 
potential  energy  of  tha  water  column.  Zt  ia  in  aaaanca  tha  same  balance  as  tha  one 
propoeed  by  Kraus  and  Turner  (1967).  In  the  sequel  tha  simplifying  assumption  will  ba  made 
that  the  affect  of  tha  diasipation  ia  reflected  in  tha  value  of  tha  mixing  efficiency 


parameter  a  (■••  also  Millar  and  Kraus,  1977).  An  expression  for  the  convactlva  term 
0  _ 

ug  /  w'T'd*  can  ba  darivad  from  Bq.  (1).  undar  tha  "slab  assuaptlon*  that  tha 
-h 

horlsontal  velocity  Is  unifora  ovar  tha  alxad  layar  depth  tha  resulting  energy  aquation  iai 
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instead  of  aaraly  prescribing  tha  net  downward  heat  flux  at  tha  surface  it  will  ba 
related  to  tha  difference  between  the  mixed  layar  temperature  and  an  apparent  atmospheric 
equilibrium  temperature  T^i 


Q  5  w'T'  (0 )  +  -  *<T  -  Ta  )  (6) 

where  TA  includes  the  effects  of  evaporation  and  solar  radiation.  The  strength  of  the 
coupling  between  the  atmosphere  and  tha  mixed  layar  is  expressed  in  the  value  of  tha  factor 
X.  Based  on  a  heat  budget  analysis  for  tonally  and  time  averaged  conditions  Raney  (1971) 
has  shown  that  it  is  of  the  order  of  8  x  10_6ms”1 . 

Imposing  a  heat  flux  in  the  form  (8)  results  in  an  extra  nonlinearity  in  the  right 
hand  side  of  the  energy  balance  (5).  However,  if  we  try  to  avoid  that  difficulty  by 
applying  just  a  given  Q  it  turns  out  that  we  throw  away  the  frontoganetic  mechanism. 

This  will  be  shown  in  the  next  section. 

The  dynamics  in  the  mixed  layer  can  for  instance  be  described  by  e  simple  Bkman 
balance.  For  the  cases  considered  here  the  mixed  layer  thicknesses  will  be  small  enough  so 
that  the  geostrophic  contribution  to  the  mixed  layer  transport  can  be  neglected.  Therefore i 
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(8) 


where  tx,t^  are  the  eastward  and  northward  components  of  the  surface  wind  stress 
respectively,  and  f  is  the  Coriolis  parameter. 


3.  Generation  of  fronta  in  the  mixed  layer 


In  the  caae  of  convergence  of  the  transport  in  the  mixed  layer,  with  associated 
downweiling,  and  under  steady  atate  conditions  Eq.  (3)  reduces  to  a  linear  balance  between 
local  atmospheric  heating  (or  cooling)  and  horizontal  advective  cooling  (heating).  The 
problem  can  then  be  solved  without  knowledge  of  the  oceanic  interior.  Also  the  mixed  layer 
depth  disappears  from  the  heat  balance  and  doesn't  need  to  be  resolved.  This  is  the  model 
studied  by  welander  (1981).  He  showed  that  if  the  mixed  layer  is  very  weakly  coupled  to 
the  atmosphere  a  discontinuity  in  the  temperature  field  is  generated  at  the  location  of 
maximum  downwalling.  The  condition  of  extremely  weak  atmosphere  -  mixed  layer  coupling 
(i.e.  X  very  small  in  (6))  is  essential  so  that  the  strong  gradients  that  result  from  the 
advection  toward  the  convergence  region  are  not  wiped  out  by  the  local  atmospheric  heating. 

De  Szoeke  (1980)  also  studied  frontoganesis  due  to  convergence  of  Ekman  transports  in 
the  mixed  layer.  His  model  is  time  dependent  and  concentrates  on  the  interplay  between 
wind  driven  entrainment  and  up-  or  downwelling.  There  is  no  surface  buoyancy  flux  so  the 
dissipative  effect  through  local  atmospheric  heating  (or  cooling)  is  fully  absent  in  his 
model.  For  variations  on  a  seasonal  time  scale  the  governing  equations  follow  from  (3) 
through  (8)  by  taking  X  •  0.  Cuahman-Roisin  (1981)  slightly  extended  De  Szoeke' s  analysis 
by  allowing  water  masses  to  escape  laterally  out  of  the  convergence  regions  (confluence). 
Also  in  that  study  the  surface  buoyancy  flux  is  neglected.  However,  it  is  well  known  that 
in  almost  all  cases  the  heat  flux  at  the  air-water  interface  is  one  of  the  dominant  factors 
that  control  the  mixed  layer  (see  for  instance  Gill  and  Niiler,  1972)  and  it  can  certainly 
not  be  neglected.  In  fact,  this  is  one  of  the  main  reasons  for  the  rather  successful  use 
of  one  dimensional  mixed  layer  models  in  simulations  of  variations  on  the  diurnal  and 
seasonal  time  scale  of  the  temperature  field  in  the  open  ocean  or  in  lakes  (Niiler  and 
Kraus  (1977),  Harleman  (1982)). 

In  the  remainder  of  this  study  the  important  role  that  can  be  played  by  the  surface 
buoyancy  flux  in  both  the  generation  and  the  dissipation  of  fronts  shall  become  obvious. 

To  insolate  the  frontogenetic  mechanism  in  its  purest  form  we  will  apply  a  wind  stress  that 
is  spatially  uniform.  Divergence  or  convergence  of  EVman  transport  is  therefore  absent. 
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Also  the  apparent  atmospheric  temperature  will  b«  uniform  in  space.  It  will  ba 

monotonically  increasing  as  a  function  of  time,  indicating  that  we  are  considering  the 
heating  part  of  the  seasonal  cycle,  when  the  seasonal  theraocline  is  formed.  Under  these 
conditions  (3)  through  (6)  simplify  to: 

2mu^ 

h  “  ag\(T  -  T)  (9) 
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where  U  and  V  can  for  instance  be  given  by  (7)  and  (8).  In  calculating  solutions  of  an 
initial  or  an  initial-boundary  value  problem  of  (10)  it  must  be  remembered  that  it  is  only 
valid  if  Sh/St  <  0.  This  condition  can  easily  be  checked  through  (9)  and,  with  >  T, 
it  is  equivalent  to: 
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Su, 

St 


u,  S(Ta  -  T) 


<  0 


It  shows  that  (in  this  model)  detrainnent  can  result  if  there  is  a  sufficient  relaxation  of 
the  wind  stress  or  if  the  surface  buoyancy  flux  increases  at  a  fast  enough  rate. 

Eq.  (10)  is  a  quasilinear  hyperbolic  partial  differential  equation  and  can  be  solved 
by  the  method  of  characteristics  (e.g.  Courant  and  Friedrichs  (1948)).  The  effect  of  the 
nonlinearity  on  the  left  hand  side  can  most  easily  be  seen  if  we  apply  it  to  a  case  with  no 
meridional  transport  (V  -  0).  Eq.  (10)  can  then  be  interpreted  as  the  ordinary 
differential  equation  (of  Ricatti-type): 


dT 

dt 
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along  characteristics  that  are  the  solution  curves  of: 
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Let  an  initial  t— pa ratura  diatribution  ba  given  by 


T(x,0)  -  Tg(x)  .  M3) 

Different  valuas  of  T0  propagata  with  different  velocities  dx/dt  aa  givan  by  (12).  Tha 
simplest  case,  with  TA  apatiaily  uni fora  and  U  and  u*  constants,  shows  that  this 
propagation  apaad  is  Xargaat  in  ragiona  with  tha  coolast  nixad  layar  water,  for  instanca 
naar  a  coastal  upwelling  region.  In  this  latter  important  application  tha  propagation  is 
directed  offshore,  with  tha  largest  spaed  naar  tha  coast,  decreasing  toward  tha  interior. 
This  results  in  a  compression  of  that  part  of  tha  initial  profile  and  under  certain 
conditions  tha  temperature  gradient  can  sharpen  so  strongly  that  a  temperature 
discontinuity  results:  a  front  has  been  generated.  It  is  this  nonlinearity  in  the 
advective  term  that  can  lead  to  such  a  steepening  of  the  gradient  with  possible 
frontogenesis.  However,  whether  a  front  will  be  generated  is  also  dependent  on  the  damping 
effect  that  is  provided  by  tha  local  heating,  represented  in  the  aquations  by  the  right 
hand  side  of  (11).  The  heating  acts  most  efficient  in  the  coolast  regions  thus  working  to 
wipe  out  an  Initial  horisontal  temperature  gradient.  These  two  processes:  local  heating 
and  differential  propagation  of  the  initial  profile  are  therefore  competitive. 

To  investigate  if  and  under  what  conditions  fronts  may  form  we  shall  solve  (11) 
through  (13)  for  a  TA  that  increases  linearly  in  time: 

Ta  -  a  +  bt  (b  >  0)  . 

The  transport  U  and  the  shear  velocity  u*  are  taken  constant.  Ondar  these  conditions 
the  solution  isi 


T(t,C) 


1  -  u<EW 


(14) 


along  the  characteristic  curve  (through  the  point  (5,0)  on  the  initial  line)  that  has  the 
equation: 
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and  5  acta  as  a  parameter . 

Tha  family  of  characteristics  given  by  (IS)  forma  an  envelope  fori 

1  •  ?k 

('  denotea  differentiation  with  respect  to  £)•  Such  an  envelope  forma  the  boundary  of  a 
region  in  tha  (x,t) -plane  where  the  charactariatica  overlap  leading  to  multivaluedneaa  of 
the  aolution.  If  an  envelope  la  formed  then  it  detarminea  the  poaition  of  polnta  at  which 
3T/3x  will  firat  become  infinite  and  a  front  forma*  It  ia  eaay  to  ahow  that  the 
condition  that  thia  happena  for  a  t  >  0  ia  equivalent  to 


ami'b  v2 


OTJ  >  X<T  <0)  -  T  )  +  (-— ) 


(18) 


expressing  the  fact  that  the  Initial  temperature  gradient  in  the  mixed  layer  muat  be  large 
enough  eo  that  the  damping  effect  by  the  local  heatinc  doesn't  prevent  the  temperature  wave 
to  "break".  Obviously  UTg  has  to  be  positive,  which  means  that  the  cooler  water  must  be 
upstream.  The  inequality  (18)  also  shows  that  having  the  largest  initial  temperture 
gradient  in  the  region  with  the  smallest  temperature  difference  with  the  atmosphere 
provides  the  most  favorable  condition  for  a  temperature  discontinuity  to  develop.  A  weak 
coupling  to  the  atmosphere  (X  small)  also  seems  favorable,  but  (17)  shows  that  the 
smaller  X  the  later  the  front  is  formed. 

Once  a  front  has  been  established  it  propagates  with  a  velocity  that  can  be  determined 
from  the  mass  conservation  law  in  integral  form.  This  is  worked  out  in  the  appendix.  The 
result  is  that 


dl 

dt 


(19) 
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where  t(t)  is  the  location  of  the  front  at  time  t  and  h2  is  the  depth  of  the  mixed 
layer  immediately  behind  the  front.  It  expresses  that,  in  this  model,  the  shallower  layer 
behind  the  front  simply  overtakes  the  deeper  mixed  layer  ahead  of  it  without  being  slowed 
down. 
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4.  Results  and  discussion 

In  the  preceding  sections  we  have  shown  that  if  condition  {18)  is  satisfied  an 
initially  smooth  horisontal  temperature  profile  in  the  nixed  layer  will  develop  into  a 
front.  The  frontogenesis  is  not  a  result  of  divergent  or  convergent  wind  or  temperature 
fields:  in  the  analysis  the  imposed  wind  stress  and  atmospheric  temperature  fields  at  the 
surface  are  both  uniform  in  the  spatial  coordinate.  The  essential  feature  is  that  the 
surface  buoyancy  flux  is  not  just  an  imposed  quantity  but  that  it  is  related  to  the 
temperature  difference  between  the  mixed  layer  and  the  atmosphere.  In  this  way  the  initial 
horisontal  temperature  gradient  in  the  mixed  layer,  combined  with  a  spatially  uniform 
atmospheric  temperature,  gives  a  nonuniform  surface  buoyancy  flux.  This  in  turn  produces  a 
horisontal  gradient  in  the  mixed  layer  depth,  and,  with  a  uniform  wind-driven  transport,  a 
compression  or  expansion  of  the  initial  tamparature  'wave*  results. 

Before  going  over  to  a  discussion  of  the  possible  structure  of  the  front  it  is 
important  to  find  out  if  the  formal  condition  (18)  can  be  satisfied  for  realistic  values  of 
the  different  physical  parameters  in  the  problem  and  if  so  whether  an  associated  front  is 
formed  within  a  reasonable  time  span  (as  given  by  (17)).  For  that  purpose  we  will 
calculate  a  few  examples.  The  first  one  is  that  of  a  near  equatorial  region,  say  at  6*S. 
The  applied  wind  stress  is  northward,  leading  to  offshore  Ekman  transport  and  associated 
coastal  upwelling.  In  the  upwelling  region  itself  the  above  formulated  model  is  certainly 
much  too  simplified  to  be  applicable.  The  focus  is  on  the  region  somewhat  farther 
offshore,  where  the  advection  of  relatively  cool  water  out  of  the  upwelling  area  still 
maintains  a  horisontal  density  gradient.  For  the  different  parameters  we  take  the  values: 

A  »  8.75  *  lO’W  (Haney,  1971);  m  ■  1.25  (see  Niiler  and  Kraus 
(1977)  tor  a  discussion  on  this  'mixing  efficisncy  parameter'); 

PQ  -  103kgm“3;  a  -  2  *  10"4#c"1i  b  -  6  *  lO^'Cs"1; 

T*  »  IcT’kgm-V2  and  TA(0)  -  Tq  *  1'C. 

Substitution  of  these  values  in  (18)  gives  the  condition  that  if  the  initial  horisontal 
temperature  gradient  satisfies 

Tj  <  -5.5  x  10"6,Cm-1 
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than  a  front  will  result.  For  instance,  if  Tg  -  -7.5  *  10_6,Cm~  tha  associated  tima 
tf  in  which  tha  temperature  discontinuity  is  ganaratad  follows  from  (17)  to  be 
ty  ■  4.2  *  106s  (about  48  days).  If  Tg  -  -10_5*Cm_1  this  tima  decreases  to  2.7  x  106s 
(about  31  days).  If  we  repeat  the  above  calculations  for  a  weaker  wind  streast 
T*  m  5  x  10“2k$F»“1s“2  tha  resulting  condition  is  that  Tg  <  -7.5  *  10“6,Cm-1 .  For 
instance,  if  Tj  “  -1 0-^*Cm-1  (17)  now  gives  ty  ■  2.5  x  io®s. 

Bast  west  temperature  gradients  of  the  above  calculated  order  of  magnitude  are  not 
unusual  in  the  near  equatorial  eastern  Pacific.  Therefore,  under  such  conditions  the 
frontogenetic  mechanism  as  studied  here  might  be  a  non-negligible  process  in  the  extremely 
complicated  dynamics  of  that  area.  Tha  larger  absolute  value  of  Tg  that  is  required  for 
a  weaker  wind  stress  reflects  the  fact  that  local  heating  is  more  dominant  under  weak 
windst  the  temperature  wave  can  only  ‘break*  if  its  initial  slope  is  steep  enough.  This 
notion  becomes  especially  dramatic  if  we  repeat  the  above  calculations  for  regions  in 
higher  latitudes.  With  the  same  parameter  values  as  given  above,  the  associated  condition 
for  Tg  becomes  unrealistic.  For  instance  at  30*S  with  a  10-1kgm-1s-2  wind  stress 
there  results  Tg  <  -27.5  x  10-6*Cm-1 .  This  is  caused  by  the  much  smaller  Ekman  transports 
in  higher  latitudes  resulting  in  a  dominant  role  of  the  local  buoyancy  flux  in  the  dynamics 
of  the  mixed  layer.  For  the  oceanic  environment  the  nonlinear  mechanism  for  steepening  of 
the  horisontal  temperature  gradient  as  proposed  in  this  study  is  therefore  of  possible 
importance  only  in  near  equatorial  regions. 

However,  this  doesn't  mean  that  the  phenomenon  doesn't  exist  in  middle  or  higher 
latitudes!  an  important  field  of  application  might  be  that  of  lakes  and  small  seas.  After 

a  period  of  strong  winds  along  the  length  axis  of  a  stratified  lake,  with  downwelling  and 

suppression  of  the  thermocline  at  the  downwind  end  and  upwelling  at  the  upwind  side  of  it, 
it  is  not  unusual  that  in  the  upper  layer  a  horisontal  temperature  gradient  results  of  the 
order  of  10-<*Cm-1 .  if  we  then  apply  a  wind  stress  of  5  x  io~2kgm~1s~2  and  take  the 
other  parameter  values  as  earlier,  a  front  develops  in  the  order  of  two  weeks,  over  such  a 

time  period  for  a  lake  of  say  40  km  length  and  with  the  initial  maximum  temperature 

gradient  close  to  the  upwind  shore,  the  jumps  in  temperature  and  mixed  layer  depth  appear 
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near  tha  middle  of  tha  lake  aftar  which  tha  front  travala  toward  tha  downwind  boundary. 

Such  a  development  is  pictured  in  Fig.  1,  showing  qualitatively  tha  same  avolution  aa  in 
tha  ocaanic  caaa,  but  with  tha  smaller  horisontal  langth  seals  of  tha  laka  leading  to  a 
shorter  tiaa  of  frontoganasis. 

Now  that  wa  have  shown  that  tha  frontogenatic  mechanism  acts  on  realistic  time  and 
langth  scales  for  lakes  and  for  tha  near  equatorial  regions  of  tha  ocean  it  becomes 
relevant  to  address  tha  problem  of  finding  tha  structure  of  such  fronts,  tthat  we  have 
constructed  in  tha  preceding  section  is  a  correct  "weak  solution*  of  (10)  with  initial 
condition  (13),  whara  for  tha  unique  determination  of  tha  jump  discontinuity  in  tha  mixed 
layer  temperature  (and  associated  depth)  wa  had  to  return  to  tha  original  mans  conservation 
law  (for  a  discussion  on  this  and  related  topics  see  Nhitham,  1974).  Tha  question  is  now 
what  physical  process  can  be  expected  to  lead  to  a  smoothing  of  that  step  structure  to  a 
continuous  profile  of  which  the  discontinuous  weak  solution  is  a  first  approximation,  one 
obvious  candidate  that  was  not  incorporated  in  the  modal  formulation  ie  horisontal 
turbulent  diffusion.  If  we  include  that  process  in  (1)  with  a  constant  eddy  diffusion 
coefficient  x  and  repeat  all  tha  manipulations  that  led  to  (10)  an  extra  term 
K(32T/3x2  ♦  32T/3y2)  must  be  added  to  the  right  hand  side  of  that  equation.  In  order  to 
get  an  idea  of  the  possible  importance  of  such  eddy  diffusion  the  aquation  will  be 
scaled.  For  that  let  nondimansional  variables  (denoted  by  primes)  be  defined  byi 

t  ■  Tr  ♦  W i  TA  -  a  ♦  0T^> 

2mu*L 

x  • Lx’’  v  • 

The  resulting  nondimansional  version  of  Xq.  (10)  with  horisontal  diffusion  added  (and 
omitting  the  primes)  1st 


♦  sgn(U)(TA  ”*>■}£■  M(TA  ”  T)2  *  * 


(20) 


where  the  parameters  w  and  c  are  defined  byi 


-13' 


O.U 


*j  a>  a> 

C  u  u  ■ 

o  a  a 

U  *j  u 
Vi  V  V 
O  "H  • 
*D  W  Q«  • 

<0  Vi  0) 


C  ^  01  1 
0  l  X 
*W  0*  V 
U  Q 
3  V  >,  » 
A  «  H 
•h  a»  .. 
u  <o  > 

U  -H 

to  0  *j  . 
♦w  *J 
•O  4J 

•  53 

a  <o 
a  at  o 
*o  o* 
<u 

^  M  P 


>,  o 

M  >  W 
H  41  -H 
*TS 

■H  *0  0) 
4J  01  H 
•H  44  0) 

C  <0  X 

•H 

c  8  § 

«  «  o 

«  X 
'm  0) 

0  •  e 

•H  4J 

C  X  «o 

$  *>  o 

i  -H  -H  . 
9>  *  *H 

0)  a 

4->  -  CL, 

w  >»  <0 

a;  -h  at 


the  initial  profiles  of  mixed  layer  temperature  and  depth 
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U  muiitm  the  relative  importance  of  local  heating  to  the  advective  heat  flux,  e  that 
of  diffusion  to  advection.  Characteristic  values  for  the  lake  example  are:  L  •  3  *  lo4mj 
|U|  •  5  *  10“1m^e”1/  u*  “  7  x  10"2ma-1;  d  “  3*C»  X  “  10m2s_1.  These  lead  to  the  values: 

U  “  0.5  and  C  •  1.1  x  10"2.  Consequently  (20)  is  a  singular  perturbation  problem,  and 
only  in  regions  with  very  strong  temperature  gradients  the  diffusion  becomes  a 
nonnegligible  part  of  the  balance.  In  this  way  the  front  as  a  discontinuity  in  the 
temperature  distribution  is  replaced  by  a  diffusive  boundary  layer  structure.  Eq.  (20) 
suggests  a  dimensional  measure  for  the  boundary  layer  width  to  be  CL,  which  is 
K2mu^/(ogld|0| ) .  Substituting  the  above  given  values  for  the  scale  factors  gives  an  order 
300m  width.  If  we  repeat  the  example  for  the  near  equatorial  ocean,  with  L  »  106m  and 
U  “  3.3m  a  (corresponding  to  a  5  *  10  kgm  s  wind  stress)  we  find  the  values: 
l>  “  2.7  (indicating  that  local  heating  is  relatively  more  important  in  those  regions)  and 
e  *  5  x  10-5  so  that  now  the  boundary  layer  is  considerably  narrower:  O(50m). 

Another  process  that  can  be  expected  to  become  active  when  the  gradient  steepens  is 
free  convection  as  a  result  of  gravitational  instability.  It  is  at  this  point  that  the 
dynamics  in  the  region  below  the  mixed  layer  becomes  relevant  also  for  the  case  of  a 
shallowing  mixed  layer.  When  the  mixed  layer  shallows  it  leaves  a  'trace*  behind  in  the 
interior.  Suppose  for  example  that  the  interior  dynamics  are  such  that  this  trace  is 
approximately  conserved.  When  the  front  passes  it  replaces  the  existing  mixed  layer  by  a 
much  shallower  and  cooler  one,  thus  creating  a  gravitationally  instable  situation.  The 
resulting  overturning  with  associated  vigorous  mixing  can  effectively  wipe  out  the 
discontinuities  and  replace  it  by  continuous  profiles.  In  this  way  the  effect  of  the 
instability  appears  to  be  horisontally  diffusive  and  a  naive  way  of  parameterizing  it  might 
be  through  an  apparent  eddy  diffusion  term  as  in  (20).  It  is  clear  that  the  'eddy 
diffusivity*  K  is  in  this  case,  among  other  things,  strongly  dependent  on  the  dynamics  in 
the  interior  in  relation  to  the  transport  characteristics  in  the  mixed  layer.  Both  in 


lakes  (see  for  example  tmberger  and  Hamblin,  1982)  and  In  the  (near)  equatorial  ocean  (e.g. 
Moore  and  Philander,  1977)  there  can  be  considerable  dynamic  activity  below  the  mixed 
layer.  Associated  values  for  x  will  therefore  differ  substantially  from  one  situation  to 
the  other. 

When  the  gradients  of  temperature  and  mixed  layer  depth  become  steeper  the  associated 
cross  frontal  pressure  gradient  also  increases.  As  a  result  a  baroclinic  jet  along  the 
front  will  develop.  Except  for  its  possible  role  in  maintaining  the  front  (Cushatan-Roisin, 
1981 )  it  is  also  interesting  to  observe  that  in  the  near  equatorial  case  such  a  jet  will  be 
equatorward  so  that  in  this  way  water  frost  the  coastal  upwelling  region,  while  being 
intensely  warmed,  is  transported  into  the  equatorial  cone. 

Of  course,  such  effects  can  only  be  studied  quantitatively  through  much  less 
simplified  models  than  the  one  we  have  presented  here.  The  purpose  of  this  study  has  been 
to  identify  through  the  simplest  possible  model  one  of  the  physical  mechanisms  that  can  be 
responsible  for  the  formation  of  fronts  in  the  upper  mixed  layer  of  lakes  or  the  near 
equatorial  ocean. 

Acknowledgements.  The  many  discussions  with  Prof.  Claes  Booth  about  the  dynamics  of  mixed 
layers  are  gratefully  acknowledged. 
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The  propagation  of  the  front. 

N*  start  with  tha  aasa  (or  buoyancy)  conservation  lav  in  tine  integrated  fens 

*2  0  t  X2 

/  /  (p(x,t)  -  p°(x,t  ))dsdx  -  /  /  Bdxdt  - 

-h(x,t)  tQ  *1 

to  0 

-  f  If  u(x  ,s,T)p(x  ,T)ds  -  /  u(x  ,*,T)p(x  ,T)ds}dt  (A1 ) 

tQ  -h<Xj,t)  -hlx^t) 

atating  that,  in  the  absence  of  vertical  advection  and  for  a  decreasing  nixed  layer  depth 
over  the  tine  interval  (tQ,t),  the  aasa  in  the  colunn  (x1(x2)  *  (-h(x,t),0)  changes  as 
a  result  of  the  total  buoyancy  flux  at  the  surface  over  that  ties  interval  (first  ten  on 
right  side)  and  the  net  inflow  of  aasa  across  tha  boundaries  in  x1  and  Xj. 

Let  the  position  of  the  front  be  at  x  ■  l(t),  and  take  x^  <  t(t)  <  i(tQ)  <  x2 
(corresponding  to  * westward"  advection >  u  <  0).  Let  the  nixed  layer  be  of  "slab  type"  and 
the  horisontal  transport  O  constant,  so  that  utx^a, t)  -  u/Mx^T)  for  -Mx^t)  <  z  <  0. 
To  derive  the  front  velocity  dt/dt  fron  (Al)  we  differentiate  first  with  respect  to  t, 
givings 


t(t) 


h(t_,t)(P(»",t)-P°(i",t0))  |£  +  /  ^<x't,(p(x,t)-p°(x,t0))+h  }dx 

*1 

t(t0) 

-  h(i+,O(pU+,t)-P°U+,t0n  /  (||  (p(x,t)-p°(x,t0))+h  -|*}dx 

x2  *2 

+  /  (|£  (p(x,t)-p°(x,t0)  )+h  -j^Jdx  -  /  Bdx-u(  p(x2,t)-p(x1  ,t)J 


(A2) 


i(t0) 


where 

h(jf,t)(p(»",t)  -  p°(t’,t  ))  -  li*  h(x,t)  ( p(x,t)  -  p°(x,t„)> 

x+i(t) 

x>Mt) 


and 
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htt  ,t)<p<i  ,t)  -  p°<<  ,tn>>  -  lia  h(x,t)(p(x,t)  -  p°(x,tn>)  . 

0  x*t<t) 
x<t(t) 


For  tg  ♦  t(tg  <  t)  w«  have 


and 


P°(i“,t0)  ♦  pU“,t> 


P°(t+,t0)  ♦  Pli\t>  , 

which  can  easily  ba  understood  froa  Fig.  Al  to  raault  froa  the  fact  that 
P°(t  »tQ )  ”  p°(t+,t(J) .  Zf  wa  lat  first  tg  ♦  t  in  (A2),  using  tha  above  relations,  and 
than  x^Xj  ♦  x,  using  tha  fact  that  h3p/3t  is  bounded  on  (x^i)  and  (i,x2>,  Eq. 
(A2)  reduces  toi 


h(i+,t)(p(l+,t)  -  pU~,t)>  D(  P(  l+,t)  -  P(l“,t>) 

which  finally  gives  tha  propagation  velocity  of  tha  front  am 


dt 

dt 


a 

h(  t*  ,t) 


(A3) 
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Through  an  analysis  of  a  more  dimensional  extension  of  the  bulk  mixed  layer 
model  of  Kraus  and  Turner  (1967)  it  is  shown  how  even  under  spatially  uniform 
atmospheric  conditions  an  initially  smooth  horizontal  temperature  gradient  in 
the  surface  minced  layer  can  still  develop  into  a  front.  Essential  for  it  is  the 
notion  that  the  net  downward  .heat  flux  at  the  air-water  interface  is  related  to 
the  difference  between  the  mixed  layer  temperature  T  and  an  apparent 
atmospheric  temperature  T  (Haney,  1971)  so  that  the  initial  gradient  in  T 
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corresponds  to  a  horizontal  variation  in  the  surface  buoyancy  flux.  As  a  result 
also  the  mixed  layer  depth  differs  from  place  to  place.  Depending  on  the 
direction  of  the  wind  driven  transport  this  produces  either  a  steepening  or  a 
flattening  of  the  initial  temperature  profile.  A  lower  bound  condition  for  the 
initial  horizontal  advective  heat  flux  is  derived  in  terms  of  the  initial 
surface  heat  flux,  the  stirring  by  the  wind  and  the  time  rate  of  change  of  the 
apparent  atmospheric  temperature.  If  that  condition  is  satisfied  a  front 
develops.  If  not,  then  frontogenesis  is  prevented  by  the  damping  effect  of  the 
local  atmospheric  heating.  It  is  shown  that  the  condition  can  be  satisfied  on 
the  scales  of  lakes  (or  small  seas)  and  in  near  equatorial  oceanic  regions. 
Mathematically  the  problem  boils  down  to  a  first  order  quasilinear  hyperbolic 
partial  differential  equation  that  is  solved  exactly  by  the  method  of 
characteristics. 


